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Abstract. We present an analysis of high-dispersion and high signal-to-noise ratio spectra of the 
DAZ white dwarf GALEX J1931H-0117. The spectra obtained with the VLT-Kueyen/UV- Visual 
Echelle Spectrograph show several well-resolved Si II spectral lines enabling a study of pressure 
effects on line profiles. We observed large Stark shifts in silicon lines in agreement with labo- 
ratory measurements. A model atmosphere analysis shows that the magnesium, silicon and iron 
abundances exceed solar abundances, while the oxygen and calcium abundances are below solar. 
Also, we compared the observed line profiles to synthetic spectra computed with variable accretion 
rates and vertical abundance distributions assuming diffusion steady-state. The inferred accretion 
rates vary from M = 2 x 10^ for calcium to 2 x 10^ gs^' for oxygen and indicate that the accre- 
tion flow is dominated by oxygen, silicon and iron while being deficient in carbon, magnesium 
and calcium. The lack of radial velocity variations between two measurement epochs suggests that 
GALEX J1931H-0117 is probably not in a close binary and that the source of the accreted material 
resides in a debris disc. 
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INTRODUCTION 

GALEX J193156.8-I-01 1745 (GALEX J193 l-i-Ol 17, thereafter) is a hydrogen-rich white 
dwarf (Vennes, Kawka, & Nemeth 2010) showing optical heavy-element lines and an 
infrared excess. The original low-resolution spectrum obtained with the New Technol- 
ogy Telescope (NTT) at La Silla Observatory showed a strong Mg llA4481 doublet and 
weaker silicon lines. Follow-up echelle spectroscopy obtained with the Very Large Tele- 
scope (VLT)-Kueyen enabled a detailed abundance study. The near-solar (±0.5 dex) 
abundances of oxygen, magnesium, silicon, calcium and iron bear the signature of an 
external supply of material accreting onto the surface of the white dwarf. Based on avail- 
able data, Vennes et al. concluded that the supply may originate from a close, sub-stellar 
companion or from a cool debris disc. 

The presence of heavy elements in hydrogen-rich white dwarfs has variously been in- 
terpreted as intrinsic to the star, or as extrinsic, i.e., supplied by the interstellar medium 
(Dupuis, Fontaine, & Wesemael 1993), by a nearby companion as in post-common en- 
velope systems (Debes 2006; Kawka et al. 2008), or by a debris disc (Zuckerman et al. 
2003; Kilic et al. 2006; Farihi, Zuckerman, & Becklin 2008). However, accretion from 
the interstellar medium is unlikely because of supply shortages (Farihi et al. 2010a). 
In the extrinsic scenarios, the elements are accreted and diffused in the atmosphere and 
envelope of the star (see Fontaine & Michaud 1979; Koester 2009). An intrinsic, or inter- 



nal, reservoir of heavy elements is also possible, but in either scenario a self-consistent 
solution of the diffusion equation must explore the effect of radiative acceleration on 
trace elements (Chayer, Fontaine, & Wesemael 1995; Chayer et al. 1995). 

As a class, the polluted DA white dwarfs, or DAZs, are often defined by the detection 
of the Call H&K doublet in optical spectra (Zuckerman et al. 2003; Koester et al. 
2005). Exceptionally, MgllA4481 is, so far, only detected in a handful of warm white 
dwarfs (but see Kawka et al. 2010) such as EG 102 (Holberg, Barstow, & Green 1997), 
GALEX J1931+0117, and two warm white dwarfs from the Sloan Digital Sky Survey 
(SDSS) that show evidence of dusty and gaseous discs (Gansicke et al. 2006; Gansicke, 
Marsh, & Southworth 2007). The presence of a large concentration of magnesium in 
the last two objects helped establish a strong link between heavy element pollution and 
dense circumstellar environments. Moreover, an infrared excess that cannot otherwise 
be explained by a cool companion, may be attributed to a dust ring as in the case of the 
DAZ white dwarf G29-38 (Graham et al. 1990). 

We present new high-dispersion spectroscopic observations that help elucidate the na- 
ture of the peculiar abundance pattern in GALEX J193 l-i-Ol 17. We revise our abundance 
measurements and explore the effect of vertical abundance inhomogeneities. 

OPTICAL SPECTROSCOPY AND ABUNDANCES 

Following the identification of GALEX J 1931-1-01 17 as a DAZ white dwarf we obtained 
a series of echelle spectra (R ^ 46000) using standard settings with the UV- Visual 
Echelle Spectrograph (UVES) attached to the VLT-Kueyen. The first series of spectra 
were obtained on UT 2009 Nov 12 (epoch 1, thereafter) with the dichroic #2 and the 
"HER_5" and "BK7_5" filters on the blue and red arms, respectively. The blue spectrum 
was centred at 4370A and covered the range 3757-4985A, and the red spectra were 
centred at 7600A and covered the ranges 5698-7532 and 7660-9464A. The first series 
of spectra were analysed by us (Vennes, Kawka, & Nemeth 2010). The second series 
of spectra were obtained on UT 2010 Mar 15 (epoch 2, thereafter) with the dichroic #1 
and the "HER_5" and "SHP700" filters on the blue and red arms, respectively. The blue 
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spectrum was centred at 3900A and covered the range 3290-45 1 8 A, and the red spectra 
were centred at 5800A and covered the ranges 4788-5750 and 5839-6808A. 

The Ha radial velocity measured at both epochs are identical within errors. We 
adopted the average vua = 37.2 ± 0.6 km s~^ The lack of radial velocity variations 
rules-out the presence of a close companion and supports the debris disc model. Figure 1 
shows Sill lines with the velocity scale centred on vna- Stark shifts (Gonzalez et 
al. 2002; Lesage 2009) are apparent in some lines. Figure 2 compares observed line 
spectra to the model line profiles that include or exclude Stark shifts. The Stark shifts 
result in a line asymmetry and a notable velocity shift. Details of the calculations are 
given by Vennes, Kawka, & Nemeth (2011). Table 1 summarizes our new abundance 
measurements based on a critical review of atomic data (Vennes, Kawka, & Nemeth 
2011). 
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FIGURE 1. Set of eleven Si II lines obtained with the VLT/UVES. The line profile analysis reveals an 
abundance of silicon nearly a factor of two above solar. The SiiiAA5041,5055, 5957,5978 lines show 
pressure (Stark) shifts in agreement with laboratory measurements. 



ACCRETION AND DIFFUSION IN THE ATMOSPHERE 



In the accretion/diffusion model the abundance at each layer is obtained by solving si- 
multaneously the mass-continuity and diffusion equations. The mass-continuity equa- 
tion for element "2" is 

M2 = 47:R^ {A2n2)vi^2 

where R is the stellar radius, and M2, A2, n2, and vi,2 are the mass accretion rate, the 
atomic weight, the number density, and the diffusion velocity of element "2", respec- 
tively. The diffusion velocity is given by the diffusion equation. Assuming diffusion 
steady-state the mass-continuity equation is solved for the abundance ratio 

C2 = — «iM2— — 2- 

where C2 << 1 is the abundance of trace element "2" relative to the main constituent 
"1" assumed to be hydrogen and Ti^2 is the diffusion time scale of element "2" related 
to the diffusion velocity by 
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FIGURE 2. Details of Sin line profiles showing an asymmetry and a velocity shift. Model including 
pressure (Stark) shifts (full Unes) are markedly superior to models excluding it (dashed Unes). 



where m and p are the total mass loading and density in the atmosphere. The selective 
radiation pressure and concentration gradients are negligible in the present case and the 
diffusion velocity is given by 



where D12 is the diffusion coefficient of element "2" which is a function of temperature, 
density and atomic charge. The diffusion coefficient for ionized species is given by 
Fontaine & Michaud (1979) and includes a correction factor for thermal diffusion. 
Montmerle & Michaud (1976) proposed to calculate the diffusion coefficient for neutral 
species based on the hard-sphere model: 



where T is the temperature in K, Hp is the proton density in cm ^, and Oq is the atomic 
radius (in A) taken as 0.91, 0.65, 1.72, 1.46, 2.23, and 1.72 for C, O, Mg, Si, Ca, and Fe, 

respectively. However, Michaud, Martel, & Ratel (1978) showed that the neutral-proton 
interaction is dominated by the induced electric dipole and that the hard-sphere model 
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TABLE 1. Abundances and accretion rates in 



GALEXJ1931+0117. 
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n(X)/n(H) 




[X/H] a 


logM(gs 1) 


[X/Si] 
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< 7 X 10-5 




< -0.54 


< 8.0 


-0.8 
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2.4 ±0.3 X 10" 


-4 


-0.28 


9.30 + 0.05 


+0.1 


Mg 


3.8±0.6x 10" 


-5 


+0.05 


7.75 + 0.06 


-0.5 


Si 


5.8±1.0x 10" 


-5 


+0.25 


8.25 + 0.07 




Ca 


7.7±0.8x 10" 


-7 


-0.42 


6.35+0.04 


-0.9 


Fe 


3.7±0.8x 10" 


-5 


+0.12 


8.50 + 0.09 


+0.0 



a [X/H]=logn(X)/n(H)-log«(X)/n(H)Q 
[X/Si]= logn(X)/n(Si)|flo„ -logn(X)/n(Si), 



may overestimate the diffusion coefficient by a large factor. The polarization diffusion 
coefficient is approximated by 



Di 2(poI) = 3.3 X IO'^tJ^^- cm^s-i 

where Ud is the dipole polarizability (in A^). We adopted averaged values of 1.68, 0.86, 
10.6, 5.49, 23.6, and 8.95 for C, O, Mg, Si, Ca, and Fe, respectively from the tabulation 
of P. Schwerdtfeger (http://ctcp.massey.ac.nz/dipole-polarizabilities). For example, at a 
reference temperature of 20 x 10^ K, the ratio of hard-sphere to polarization diffusion 
coefficients is £)i,2(h — s) /Di 2(pol) ^ 3 for Mg, Si, Ca, and Fe, ^ 5 for C, and 6 for 
O. The diffusion coefficient averaged over ionization species is given by: 

i 

where xi = n2j/n2 are the ionization fractions for element "2" and ^1,2,/ are the corre- 
sponding diffusion coefficients. 

Figure 3 shows the predicted abundance profiles of O, Mg, Si, and Ca in the atmo- 
sphere of GALEX J 193 1+0 117 calculated using the procedure described above. The 
accretion rates were set at the best-fitting values. The effect of neutral species on the 
profiles of Mg, Si, and Ca is negligible because of their relatively high ionization frac- 
tions. On the other hand, the predicted oxygen abundance profile in the line forming 
region is dominated by diffusion of neutral oxygen atoms. Employing the hard- sphere 
diffusion coefficient for neutral oxygen leads to an overestimation of the diffusion ve- 
locity and an underestimation of the abundance by up to a factor of 4 relative to more 
realistic calculations that employ the polarization diffusion coefficient as determined by 
Michaud, Martel, & Ratel (1978). 

Table 1 lists the accretion rates M for each element measured by fitting spectral 
syntheses based on calculated abundance profiles to the observed line spectra. The 
inferred abundances in the accretion flow reveal a relatively high abundance of oxygen, 
silicon, and iron relative to carbon, magnesium and calcium. 
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FIGURE 3. Abundance profiles in tlie line forming region predicted by the accretion-diffusion model. 
The accretion rates tiiat fit the observed O, Mg, Si, and Ca line spectra are given in g s~ ^ The abundances 
measured using homogeneous model atmospheres are shown with dotted horizontal lines. Diffusion pro- 
files excluding neutral species (short dashed lines) overestimate the abundances in the line forming region 
relative to models including neutral species and realistic polarization diffusion coefficients (full lines), 
while models using hard-sphere diffusion coefficients (long dashed fines) underestimate the abundances. 
Mass densities corresponding to optical depth uruty at A = 7400 (leftmost), 4700, and 1260 A (rightmost) 
are marked with vertical arrows. 



SPECTRAL ENERGY DISTRIBUTION 

A lack of radial velocity variations suggests that the source of the accreted mate- 
rial is a circumstellar debris disc. Simulated infrared spectral energy distributions for 
GALEX J193 1-1-01 17 including a spectral synthesis of the white dwarf added to a warm 
disc (Fig. 4) show that the observed distribution (Vennes, Kawka, & Nemeth 2010) can 
be matched by a warm disc (900 K) or a hot disc (1650 K) near sublimation temperature. 
The inferred dust temperature must be below sublimation temperatures (see Lodders 
2003) for the main constituents such as calcium (1659 K), or silicon and magnesium 
(1354 K). The corresponding disc sizes are 90Q.wd for the hot disc, and, compensating 
for its reduced emissivity per unit area, 2200D.WD for the warm disc, where Q.wD is the 
solid angle subtended by the white dwarf itself. The observed spectral energy distribu- 
tion extends to the 2MASS measurements in the infrared and measurements at longer 
wavelengths are required to further constrain the model. 

A comparison with other warm DAZ white dwarfs supports these conjectures. The 
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FIGURE 4. Predicted spectral energy distribution of GALEX J193 1+01 17 (full lines) extending beyond 
the JHK photometry and including the emergent flux of a hot white dwarf (DAZ) added to black-body 
emission from a disc at 900 or 1650 K. The white dwarf dominates the ultraviolet region (GALEX N^y 
and Fuy) while the disc is predicted to dominate the infrared region. 



DAZ PG1015+161 (Koester & Wilken 2006) is similar to GALEX J1931+0117, and, 
based on JHK and Spitzer observations. Jura, Farihi, & Zuckerman (2007) inferred the 
presence of a disc with an inner temperature ranging from 800 to 1000 K and emitting 
area from 700 to 1200 Q.wd- In a similar case, Farihi, Jura, & Zuckerman (2009) inferred 
the presence of a disc with a blackbody temperature of 1500K around PG1457-086. The 
great diversity in disc emissivity is illustrated by the case of HS0047+1903. Despite 
having a similar calcium abundance, no evidence of a disc is found in JHK and Spitzer 
observations (Farihi et al. 2010b). 

CONCLUSIONS 

We presented a model atmosphere analysis of the high-metallicity white dwarf 
GALEX J1931+0117. The abundance pattern obtained using homogeneous model 
atmospheres shows that magnesium, silicon, and iron are near or above solar abun- 
dances, while carbon, oxygen and calcium are below solar abundances. However, a line 
profile analysis performed using vertical abundance distributions obtained by solving 
the steady-state diffusion equation shows that the accretion flow is rich in oxygen, 
silicon, and iron while it is depleted in other elements. Although the oxygen abundance 
is below solar in the line forming region, it must be supplied in larger quantity because 



of its short diffusion time-scale relative to other elements. 

The effect of pressure shifts are apparent in several strong silicon lines. This effect 
predicted by Hammond (1989) and Krsljanin & Dimitrijevic (1992) also impacts radial 
velocity measurements. A lack of radial velocity variations between two epochs that are 
123 days apart rules out a close binary scenario for the origin of the accreted material. 
We are left with the possibility that the infrared excess belongs to a warm disc of debris 
material that accretes onto the white dwarf surface. 

In conclusion, our analysis of the abundance pattern in GALEX J1931+01 17 and the 
absence of a close companion, as well as a comparison with similar cases (e.g., DAZ 
PG 10 15+ 161) support the likely presence of a dusty disc around GALEX J 193 1+01 17. 
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